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Abstract 

We study the radial- hedgehog solution in a three-dimensional spherical droplet, with homeotropic 
boundary conditions, within the Landau-de Gennes theory for nematic liquid crystals. The radial- 
hedgehog solution is a candidate for a globally stable configuration in this model framework and is also 
a prototype configuration for studying isolated point defects in condensed matter physics. The static 
properties of the radial-hedgehog solution are governed by a nonlinear singular ordinary differential 
equation. We consider two different limits separately - the vanishing core limit and low-temperature limit 
respectively. We use a combination of Ginzburg-Landau techniques, perturbation methods and stability 
analysis to study the qualitative properties of the radial-hedgehog solution, both in the vicinity of and 
away from the defect core. We establish the instability of the radial-hedgehog solution with respect to 
biaxial perturbations in certain parameter regimes and demonstrate the stability of the radial-hedgehog 
solution in other parameter regimes. Our results complement previous work in the field, are rigorous 
in nature and give information about the role of geometry and temperature on the properties of the 
radial-hedgehog solution and the associated biaxial instabilities. 



1 Introduction 

Defect structures have attracted a lot of interest in the liquid crystal community [20, 22, 23] [21]. Defect struc- 
tures in liquid crystalline systems are usually modelled within the Landau-de Gennes framework, whereby 
the liquid crystal configuration is mathematically described by a symmetric, traceless 3x3 matrix, known as 
the Q-tensor order parameter [5] . The Q-tensor can be written in terms of its eigenvalues and eigenvectors 
as shown below - 
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Q = ^A i e i (g)e. i , ^ X t = (1) 
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where Xi are the eigenvalues and are the corresponding orthonormal eigenvectors. The liquid crystal is 
said to be in the (i) isotropic state when Ai = for i = 1 . . . 3, (ii) uniaxial state when Q has a pair of equal 
non-zero eigenvalues and (iii) biaxial state when Q has three distinct non-zero eigenvalues [19j . 

A prototype example of such a confined system is a spherical droplet with strong radial anchoring or 
homeotropic (normal) boundary conditions. This example has been widely studied in the literature, especially 
from a numerical point of view, and it is generally believed that there are two competing equilibrium 
configurations - (a)the radial-hedgehog solution which has a single isolated point defect at the droplet centre 
and (b) the biaxial-torus solution where the point defect broadens out to a ring-like structure around the 
droplet centre [9j [T3l [TTl [23j [24]. The radial- hedgehog solution is purely uniaxial everywhere except for 
an isotropic point at the droplet centre whereas the biaxial-torus configuration exhibits a high degree of 
biaxiality around the droplet centre. The isotropic point in the radial-hedgehog solution and the biaxial 
ring in the torus solution are interpreted as being defect structures since they are localised regions of abrupt 
changes in the eigenvalue structure. 



This paper aims to build a self-contained mathematical description of the radial-hedgehog solution within 
the Landau-de Gennes framework. Firstly, this is an interesting mathematical problem in its own right since 
the radial- hedgehog solution is a rare example of an explicit solution of the Landau-de Gennes Euler-Lagrangc 
equations in (|T5|). Moreover, the corresponding scalar order parameter is a solution of an ordinary differential 
equation (see (l20l) 'l and hence, has a tractable and yet non-trivial mathematical structure. Indeed, this is 
the first step in the mathematical theory of defects in liquid crystalline systems. Secondly, a systematic 
mathematical analysis of the radial- hedgehog solution is crucial for understanding the structure and locations 
of point defects in liquid crystalline systems, the multiplicity of uniaxial solutions and the characterization 
of the competing biaxial structures. 

To further elaborate on the above, radial-hedgehog solutions can be thought of as prototypical vortices 
in the Ginzburg-Landau theory for superconductors 2 . More precisely, the radial-hedgehog solution can 
be interpreted as being a degree +1 vortex in three dimensions. There is a very well-developed theory for 
the structure, location, multiplicity and stability of vortices in Ginzburg-Landau theory, especially in two 
dimensions but generalizations to higher dimensions are non-trivial HH1 [HI UH] • One of the main aims 
of this paper is to clearly demonstrate the analogies and differences between the mathematical formulation 
of radial-hedgehog solutions in the Landau-de Gennes framework and Ginzburg-Landau vortices. Once the 
inter-relationship is correctly understood, this will contribute to a sound theoretical foundation for defects 
in liquid crystals and Ginzburg-Landau numerical methods can also be used for the simulation of defects in 
liquid crystalline systems. We deal with two separate limits in this paper - the low-temperature limit where 
the governing ordinary differential equation has an almost Ginzburg-Landau structure and the vanishing 
core limit where there are important technical differences. In particular, we cannot exploit Ginzburg-Landau 
techniques to describe the isotropic defect core in the vanishing core limit. The low-temperature limit is 
relevant for liquid crystalline systems deep in the nematic phase where we expect to see a high degree of 
orientational ordering. The vanishing core limit is relevant for materials whose elastic constants are typically 
much smaller in magnitude than the thermotropic parameters and quoted values in the literature suggest 
that this limit is relevant for commonly used liquid crystalline materials [21]. More generally, although the 
study of uniaxial states can be viewed as a generalized Ginzburg-Landau theory from R 3 — > R 3 [16 , biaxiality 
presents a whole host of new mathematical challenges, outside the scope of Ginzburg-Landau theory |14) . In 
particular, there is no analogue of a biaxial instability in the current Ginzburg-Landau literature and such 
instabilities play a pivotal role in Landau-de Gennes theory. 

The paper is organized as follows. In Section [21 we prove the existence of a radial-hedgehog solution in 
the Landau-de Gennes framework and establish bounds for the corresponding scalar order parameter. In 
Section[31 we derive a series expansion for the radial-hedgehog solution near its isotropic core and demonstrate 
its similarity with three-dimensional vortices in Ginzburg-Landau theory [5] . We then show that the radial- 
hedgehog solution cannot be a global Landau-de Gennes energy minimizer for sufficiently large droplets, 
for sufficiently low temperatures by means of an explicit comparison argument . This result is qualitatively 
similar to a result reported in [9] but our method of proof is different. In Section [4j we focus on the 
low-temperature limit and the resulting Ginzburg-Landau structure of the governing ordinary differential 
equation. We use shooting arguments to establish qualitative properties of the corresponding scalar order 
parameter and use Ginzburg-Landau techniques to prove the uniqueness of the radial-hedgehog solution in 
this limit. We deal with the vanishing core limit in Section [SJ whereby the governing ordinary differential 
equation does not have a Ginzburg-Landau structure. Our results are weaker in this case and describe the 
far-field properties, away from the defect core. In Section [6[ we perform a linear stability analysis of the 
radial-hedgehog solution and this stability analysis gives insight into the effect of the ball radius on the 
associated equilibrium structure. In Section [71 we discuss our results and how they complement previous 
work in this area. 
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2 Preliminaries 



We study the qualitative properties of radial-hedgehog solutions inside spherical droplets, B(0,R) C R 3 , 
where 

B(0,i?) = {reIR 3 ;|r|<i?} (2) 

and R > is independent of any model parameters, subject to strong radial anchoring conditions. We work 
within the Landau-de Gennes theory for nematic liquid crystals, in the low-temperature regime. 
Let S C M 3x3 denote the space of symmetric, traceless 3x3 matrices i.e. 

§ d = f {Q G M 3x3 ; Q l3 = Qji, Q u = 0} 

where we have used the Einstein summation convention; the Einstein convention will be used in the rest of 
the paper. The corresponding matrix norm is defined to be 

|Q| = vW= VQ^j '../ L...3. 

We recall from [UJ [TO] that an arbitrary Q G S can be written as 

Q = s ( n ® n 1 ] + r ( m ® m 1 



3 / V 3 

where n, m are orthonormal eigenvectors of Q, s, r are scalar order parameters and we either have < r < | 
or § < r < 0. If Q G S is uniaxial, then this representation formula can be simplified to 

Q = s (g) n — -I 

where n is the leading eigenvector of Q and s is a scalar order parameter that measures the degree of 
orientational ordering about n. 

The Landau-de Gennes energy functional is given by [H] 

ZMQ] - / §|VQ| 2 + / S (Q) dV (3) 

Jb(o,r) 1 

where |VQ| 2 = Ylijk=i (l^f) ^ s ^ ne elastic energy density, L is a small material-dependent elastic 
constant and fs is the bulk energy density given by 

MQ) = -ytrQ 2 - |W + t (trQ 2 ) 2 + C (a 2 , 6 2 ,c 2 ) . (4) 

The form ([4]) is the simplest form of the bulk energy density that allows for a first-order nematic-isotropic 
phase transition; here b 2 , c 2 are material-dependent positive constants and a 2 > is a temperature-dependent 
parameter. Typical values of these characteristic constants are a 2 = 0.042 x 10 6 (T* — T)N/ Km 2 ,b 2 = 
0.64 x 10 6 7V/m 2 ,c 2 = 0.35 x 10 6 iV/m 2 where T is the absolute temperature and T* is a characteristic 
temperature below which the isotropic phase Q = ceases to be a locally stable stationary point of /b in (HJ) 
[191 [21] . We work in the low-temperature regime where the bulk energy density attains its global minimum 
on the set of uniaxial Q-tensors given by [TO] 

Qmin = |Q e S, Q = s+ fnigin- ^IJ 1 (5) 
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with n <E S 2 and 



b 2 + V& 4 + 24a 2 c 2 
S+ = Ac 2 ' (6) 

In particular, as a 2 increases, we move to lower temperatures deep in the nematic phase. The additive 
constant C (a 2 , b 2 , c 2 ) in @ ensures that /s(Q) > for all QeS. The admissible space is defined to be 

Aq = {Q G VK 1 ' 2 (5(0, 5); 5) ; Q = Q 6 on 95(0, R)} (7) 
and the Dirichlet boundary condition Q(, e Q m in is specified to be 

Q i = * + (i®i-|l). (8) 

This is referred to as strong radial anchoring in the liquid crystal literature |17l [§] , since A is the unit vector 
in the radial direction. The physically observable, equilibrium configurations correspond to either global 
or local minimizers of leg in Aq. For completeness, we recall that W 1 ' 2 (5(0, R); S) is the Sobolev space 
of square-integrable Q-tensors with square-integrable first derivatives [7]. The corresponding M^ 1,2 -norm is 

/ \ 1/2 

given by ||Q|j w 1 < 2 (B(o,H)) — \Sb!qr) IQI 2 + I^QI 2 dx\ . In addition to the W 1 ' -norm, we also use the 

L°°-norm in this paper, defined to be ||Q||ioo(B(o,fi)) = ess su PxeB(o.fl)IQ( x )l- 

In what follows, we consider two different limits: the L — > limit which is referred to as the vanishing 
core limit and the a 2 — > oo limit, which is referred to as the low-temperature limit. The reason for making 
this distinction will become clear in the subsequent sections. We work in a dimensionless framework and as 
outlined in |171 [§] , we introduce the following dimensionless variables : - 



where £ = \j 21 fi L ■ It is straightforward to show that the corresponding dimensionless energy density is 

e(Q, VQ) = i|VQ| 2 - ^trQ 2 - V6trQ 3 + ~ (trQ 2 ) 2 + C[t) (10) 

where t = > is a dimensionless reduced temperature, C(t) is an additive constant that ensures 

hit) = -ftrQ 2 - VfStrQ 3 + \ (trQ 2 ) 2 + C(t) > and t > 1 throughout the paper. 
We employ a second change of variable 

Q = ^; r = Vtf (11) 
h+ 

where 



K = !±^±5. (12) 

This choice of the dimensionless variables is especially relevant for the low-temperature limit, as will be 
demonstrated in Section 2] The corresponding dimensionless energy density is 

e(Q, VQ) = i|VQ| 2 - itrQ 2 - ^+trQ 3 + |- (trQ 2 )' . (13) 
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One can readily compute the Euler-Lagrange equations associated with the energy functional, 

Zlg[Q] = / e(Q,VQ) dV, (14) 
Jb(o,r) 

where R = y^|. In what follows, we drop the tilde on the dimensionless variables for brevity and all 
subsequent results are to be understood in terms of the dimensionless variables. The associated Euler- 
Lagrange equations are [16j E] - 



AQij = -Qy ( QifeQfcj - 3 tr (Q )J + -pQiitr(Q^), i,j = 1,2,3, (15) 

where the term ^-tr(Q 2 ) is a Lagrange multiplier associated with the tracelessness constraint. It follows 
from standard arguments in elliptic regularity that any solution Q of the nonlinear elliptic system (1151) is 
smooth and real analytic on 5(0, i?) |14j . In particular, all global and local energy minimizers in Aq are 
classical solutions of (fT5|) . 

Radial-hedgehog solutions are examples of spherically-symmetric uniaxial solutions of the system (1151) in 
the admissible space .Aq and have the form 



Here the scalar order parameter h only depends on the radial distance r — |r| from the origin and the 
corresponding admissible space is defined to be 

A h = {h€ W 1 ' 2 ([0, R],R) ; h(R) = 1} . (17) 

We note that Q € W 1 ' 2 (B(0, R); S) necessarily implies that h € W 1 ' 2 ([0, R]; K) since the eigenvalues of a 
symmetric matrix are Lipschitz functions of the matrix components [25) and hence, Ah is a natural choice 
for the admissible space. There may be multiple spherically-symmetric solutions of (|15[) but we define a 
radial-hedgehog solution to be an energy-minimizing spherically-symmetric solution in the sense described 
below. 

Proposition 2.1. (a) Consider the energy functional 

m =[\*(l(^)\^ +m ) dr (18 , 



2 V dr 



r 2 



defined for functions h € Ah, where 



/(M = -y-^ 3 + §^ 4 + CW, (19) 

f(h) > for h € Ah and /(l) = 0. There exists a global minimizer h* € Ah for I in A18\) . The function h* 
is a solution of the following singular nonlinear ordinary differential equation 

d 2 h 2 dh 6h , , o 3h + o , 9N 

— + -- =- h + h 3 + -±(h s -h 2 ) (20) 

ar -2 r dr r z t 

subject to the boundary conditions 

h(0) = and h(R) = 1. (21) 
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Moreover, h* is analytic for all r > 0. 

(b) Define the radial-hedgehog solution by 

where h* is a global minimizer of I[h] in U8\) . in the admissible space Ah ■ Then Q* is a solution of the 
Landau-de Gennes Euler- Lagrange equations h!5\) i.e. is a stationary point of the Landau-de Gennes energy 
functional. Moreover, these solutions satisfy the following energy bound 



Ilg[Q*] < 127ri? (23) 

where I has been defined in {1$ . 

(d) The function h* satisfies the following bounds for r £ [0, R] - 

0<h*(r)<l re[0,R}. (24) 

Proof 

(a) Consider the energy functional I[h] denned for h £ Ah- Firstly, we note that the admissible space Ah 
is non-empty. Indeed, the constant function h(r) = 1 for r £ [0, R] belongs to Ah- Secondly, the functional / 
in (|18[) is bounded from below i.e. I[h] > for h £ Ah and is weakly lower semicontinuous on our admissible 
space (since the integrand is convex in dh/dr). The existence of a global minimizer h* £ Ah now follows 
from the direct methods in the calculus of variations [7J. 

It is straightforward to compute the Euler-Lagrange equations associated with the functional I in (fP8]l 

i.e. 

d_ ( de(h,h') \ _ de(h,h') 
~cfr I dh' J ~ dh 

where ti = dh/dr, e(h, h!) = r 2 U (^f + ^ - ^ - '-^-h 3 + ^-/i 4 ) . One can check that the correspond- 
ing Euler-Lagrange equation is indeed the ordinary differential equation in (|20p and a global minimizer h* 
is necessarily a solution of these Euler-Lagrange equations. 

The boundary condition h* (R) = 1 follows from our definition of the admissible space Ah ■ All functions 
h € Ah are necessarily continuous since h e W 1 ' 2 ([0, R], R) => h € C°' a ([0, R], R) for some < a < 1/2 
from the Sobolev embedding theorem [7]. The boundary condition h(0) = follows from the continuity of 
h*(r) for r G [0,i?]. We proceed by contradiction and assume that > kg for r £ [0, ro], for some fixed 

ho > and < ro << 1. Since h* is continuous, we deduce that h* has a fixed sign near the origin and we 
further assume that h*(r) > ho > for r £ [0, ro]. Consider the governing equation (|20| ; it can be re- written 
as 

l(r>£)=Sk + r>(k°-k + 3 ±(k>-k>)) (25) 
where h + has been defined in (fl"2"|) . Then, we have 



2<ih > I 6h(r')dr' + Cr 3 + e 2 h{e) for r £ (0,r ) (26) 



dr 

where < e < r/10 is fixed, h (e) = ^\ r=<L and C is a constant independent of e. We note that h (e) can 
be bounded independently of e i.e. |^| < C(t) for r £ [0, R] from [Tl]. Squaring both sides of (|2"6")l and 
integrating from e to r, we obtain 

[ dri>J^dt + C"r 3 + e 2 h'(e)£±dt for r £ (0,ro), (27) 
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where C is a constant independent of e. In the limit e — > 0, ([2"?} contradicts the hypothesis that h 6 
W 1 ' 2 ([0, R]; R) from which we must have 

dh\ , 
— dr < oo. 



dr J 

Therefore, we deduce that h(Q) = for any solution of (|2U1) in .A^ and /i* € ^ is a solution of ([20]) . subject 
to the boundary conditions (|2ip . The analyticity of /i* now follows from standard arguments in the theory 
of ordinary differential equations [TT] . 

(b) Given h* , define a radial- hedgehog solution as follows 

Prom [25], Q* € W 1 ' 2 (B(0, R); S) necessarily implies h* G M /1 ' 2 ([0, R]; R + ) and the preceding arguments 
necessarily require h*(0) = 0. One can directly check that 

l cg [Q*} = 4nl[h*} (28) 

and that Q* is a solution of the Euler-Lagrange equations (TT51) . since h* is a solution of the ordinary 
differential equation (|2U|) . subject to the boundary conditions (j2T1) . 

The function /i* has been defined to be the global minimizer of the functional / in (|18p . in the admissible 
space Ah- However, the constant function, h(r) = 1 for r € [0, R], belongs to Ah and hence 

I[h*] < I[h] = 3R. (29) 



The energy bound on Xcg[Q*}, where Q* = J |/i*(r) - follows from 

(c) The upper bound |/i*(?')| < 1 follows directly from a result in [15] where we establish that every 
solution Q of the system (fTS)) in the admissible space „4q satisfies the global upper bound 

IQWI < 1. 

The radial-hedgehog solution Q* is a solution of the system (|T5|) and 

|Q*(r)| = \h*(r)\ 

where r — |r|. The upper bound |/i*(r)| < 1 follows immediately. 

The lower bound h*(r) > follows from the energy minimality condition. We assume that there exists 
an interior measurable subset 

Cl = {r e Q; h*{r) < 0} C 5(0, R) 

with h*{r) — on dfl. We note that £1 must be an interior subset because of the boundary condition Qf, in 
©. We define the perturbation 

V ' \-h*(r), re!l. 1 ' 

One can then easily check that 

I\h*\ - I[h*] = f (h*f + ^±h* 3 dV = [ ^h* 3 dV < (31) 
.la 3 o Ja t 
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since h*(r) < on f2 by assumption. The inequality pip contradicts the global minimality of h* in Ah and 
hence, we deduce that h*(r) > for r € [0, i?]. The inequalities (|24p now follow. □ 

In summary, in Proposition 12.11 we prove the existence of a radial- hedgehog solution of the form l]22p. 
that can be interpreted as being a Landau-de Gennes energy minimizer within the class of radially-symmetric 
configurations. This radial- hedgehog solution satisfies the energy bound (|23]) and the corresponding scalar 
order parameter h* is bounded from both above and below as shown in (|24p . The radial- hedgehog solution 
has a single isolated isotropic point at the origin where h* vanishes and this isolated isotropic point is 
interepreted as being a defect point, since the radial-hedgehog solution is strictly uniaxial everywhere else. 
In the next section, we study the isotropic core of the radial-hedgehog solution and the manifestation of 
biaxial instabilities within this core. 



3 The Isotropic Core and Biaxial Instabilities 

Proposition 3.1. Let h* be a global minimizer of the energy functional I in H8\) . Then h* is a solution 
of the ordinary differential equation \20\) subject to the boundary conditions V21\) . As r — > 0, we have the 
following series expansion for h* 



h* (r) = V" a n r n = a 2 r 



1 h o(r 2 

14 V 



as r — > 



(32) 



where a n = for all n odd and a 2 > is an arbitrary constant. In addition, as the domain size R — > oo, we 
also have that [17, 

(33) 



2 2 

' <h*(r)<- '" 



14 



t\ 2 



2 1 



< | < 3 since t > 1. Therefore, for R sufficiently large, we have the following 
bounds on the constant a 2 in 



where X 2 t = ■ 9+8t+3VSTm ^ t 



1 



1 



1 



(34) 



Comment: Equation h3 6 2\l is identical to the series expansion for three-dimensional vortices near the 
origin, within the Ginzburg- Landau theory for superconductivity 

Comment: The limit R —¥ oo corresponds to either the L —} limit in |2P or the t — > oo limit in 113\) or 
the doubly infinite limit L — ¥ 0, t — > oo. 

Proof: From Proposition [2Tl we have that h* is analytic for r > 0. We seek a power series expansion of 
h* around the origin with h*(0) = 0, of the form 



h*(r) = J2 a nr n 0<r<i? o 



(35) 



71=1 



where Rq is the radius of convergence. 

We substitute the ansatz (|55|) into the ordinary differential equation (|2"0]) and equate the coefficients of 
r" on both sides of (|2U)) . Straightforward computations show that 



o-i — 03 = 0, a 2 > is arbitrary 

__ 

14 



«4 

h*(r) =a 2 



r 2 - 



14 



(36) 
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where 02 > since h* is non- negative from Proposition 

Next, we show that the formal expansion (|35l) involves no odd powers of r. Direct computations show 
that a\ — a 3 = 0, as stated in We proceed by induction. Suppose that a2 n +i = for n = . . .p. We 

show that a2 P +3 = too. Consider the left-hand side of the ordinary differential equation QUI) i.e. 



2u* 



d 2 h 

dr 2 



2 dh* 6h* ^ „ n r 2 -, 
--^-^T=H r 2 a n [n 2 + n-6] 



p=0 



so that the coefficient of r 2p+1 is (4p + 2) (p + 3)02^+3. We compute the coefficient of r 2p+1 on the right-hand 
side of (|2"D|) . One can directly show that 



3h± 



h* 3 - h* 2 



where 



?2p+l 



3/l, 



3(a 1 a2 P -i + ala 2p -3 + ■ ■ ■ + + a j+i 



6/1-1 



(aia 2p + a 2 a2 P _i + . . . + a p a p+1 ) 
(37) 



where the term involving 2 P +i comes into play if p £ is a positive integer. One can check (|37l) by noting 

that the coefficient of r 2p+1 in the series h* 2 is X)n=i 2a„a2 P +i_ n so that both {rt, 2p + 1 — n} < 2p + 1 
and one of {n, 2p + 1 — n} is odd. Similarly, we note that the coefficient of r 2p+1 in the series h is 
O ap+i + 5Z«=i 3a^a2 P +i-2n, {«, 2p + 1 — 2n} < 2p + 1, 2^±1 and 2p + 1 — 2n are necessarily odd. However, 

3 

from the hypothesis, a 2n +i = for n = . . .p. Therefore, &2p+i = in ([3"7]) and since 



5 2p+l 



(4p + 2)(p + 3)a 2 p + 3, 



we deduce that ei2p+3 — as required. 

The bounds (f3"3"|) have been established in [T7J[H] and are valid in the R — > oo limit. The inequalities AM]) 
follow from (|3"3"| and the limit 

h*(r) 



0,2 



lim 

i — >o r z 



□ 



Proposition 3.2. Consider the radial-hedgehog solution 

Q*(r) = J|ft*(r) 



3' 



where h* is a global minimizer of I in \18\) in the admissible space Ah- Then Q* is not the global minimizer 
of Ilq in (TT^P in the admissible space 



•Aq = < Q e W ' (B(0, R);S) ; Q 



r r 

- ® - 

r r 



on OB(0, R) 



(38) 



(which is the admissible space defined in ^ in terms of the dimensionless variables in (OJ) and for 
sufficiently large values of the domain size R and the reduced temperature t. In particular, the biaxial state 



Q(r) 



Q*(r) cr < r < R, 



0<r < a 



(39) 
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where z = (0, 0, 1) is the unit-vector in the z-direction, has lower free energy than the radial-hedgehog solution 
for 

a = 10 t = 200. (40) 



Proof: Consider a general biaxial perturbation (f3"!)| 

r Q*(r) +p(r) (z®z - §" 



Q(r) 



< r < cr 



Q*(r) cr < r < R, 



where p : [0, R] — > K is non-zero for < r < a and p(r) = for all u < r < i?,, Q* is the radial-hedgehog 
solution in ([22]). r = (x,y,z) is the position vector, z = (0,0, 1) is the unit-vector in the z-direction and I 
is the 3x3 identity matrix. In particular, the perturbation Q is localised in a ball of radius a around the 
origin or equivalently, is localised around the isotropic core of the radial-hedgehog solution and the radius a 
will be determined as part of the problem. 

Let (r,0,<j>) with r G [0, R], 9 € [0, 7r] , <f> € [0, 2tt) denote a spherical coordinate system centered at the 
origin. Straightforward computations show that 



|VQ| 2 = |VQ 



* |2 



V6 



dh* dp 
dr dr 



trQ 2 = trQ* 



3 

, 2 



p 2 (r) + V6h(r)p(r) I cos 



cos 2 9 

9- 1 - 
3 



9 X 



trQ J = trQ" + -p%r) + W - + —= \ h*{r)p z {r) cos 2 9-~)+ -h*\r)p(r) cos 2 9 - - 



V6 



1 



1 



(trQ 2 )" = (tiQ* 2 y + ^p\r)+6(h*{r)) 2 p 2 (r) [ cos J 



(h*(r)) 2 p 2 {r) + 2V6h*(r)p{r) 



Noting that 



cos o — 



cos 



h* 



sin 



and 



we obtain the following 



3 Sm 



= 



45' 



' % G [Q] -X LG [Q*] 



4tt 



- Vp 2 M"2V6^rV(r) 



r 2 h\ 
2t 



4 4/ \ i 28 2 2/ x 

gP (r) + — /i p (r, 



where h + has been defined in (jT2"j) . Recalling the bounds (f55j> . we have that 



1 

4tt 

< 



Ilg[Q] -Ilg [Q1 



< 



L. p2(r) _ 2 V6^trV(r) + ^ 



4 4/ x , 28 
9 P (r) + 15 



dr 



(41) 



(42) 



p 2 {r) } dr (43) 



10 



Where X " = 9 + 8t+3V9T8i - LCt 



One can then directly substitute (|4"4")l . cr = 10 and £ = 200 into (fT2"|) and (|4U1) to find that the associated free 
energy difference 



1 

47T 



lLG[Q}-iL G [Q* 



< o 



i.e. we have found a biaxial perturbation localised in a ball -B(0,cr), that has lower free energy than the 
radial-hedgehog solution for t — 200. Therefore, the radial- hedgehog solution cannot be a global Landau-de 
Gennes minimizer in this regime. Proposition 13.11 now follows. □ 

The instability of the radial-hedgehog solution with respect to biaxial perturbations has been theoretically 
demonstrated in [3], in the limit t — > oo and R — > oo. The authors in [5] consider the second variation of 
the Landau-de Gennes energy functional and treat the instability condition as a Schrodinger eigenvalue 
problem, which has to be solved numerically. We construct an explicit biaxial perturbation, localized near 
the isotropic core of the radial-hedgehog solution and show that this biaxial perturbation has lower free 
energy than the radial-hedgehog solution, for sufficiently low temperatures. The biaxial perturbation is 
energetically preferable only when localized in a ball -8(0, cr) centered at the origin and one can check that 

Xlg[Q] — ZlgIQ*] > if is too small or too large i.e. a needs to be large enough for the biaxiality 
to manifest itself and yet be small enough so as not to perturb the far-field properties. We demonstrate 
instability for t — 200 and a = 10 and the parameter regimes can be investigated more systematically. 
In particular, our approach in Proposition 13.11 gives insight into how to quantify the instability regime 
analytically. 



4 The limit t — ► oo 

Consider the ordinary differential equation in (120 

dr 2 

in the limit t — > oo. In the limit t - 



2 dh 6h 

' Zl 2 = ~ k 

r dr r z 
oo, 



3h, 



h±<jL 

t ~ y/i 



for some (3 > independent of t and hence for any non-negative solution h, we have 



3h H 



(h 3 - h 2 ) 



« h 



since < h(r) < 1. Recall that the upper bound in (j2~i)) applies to all solutions of (|2T)|) and not just the 
radial-hedgehog solution. In the limit t — ^ oo, the ordinary differential equations (|20[) approximately reduces 
to 

d 2 h 2dh 6h , ,„ 

(45) 



dr 2 



r dr 

3h + 



-h + h 6 



although the influence of the perturbation term 22± (/j 3 — ft 2 ) needs to be carefully quantified. The limiting 
problem (|45p has a very similar structure to the governing ordinary differential equation for vortex solutions 
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in the Ginzburg-Landau theory of superconductivity [I] . Vortex-solutions have been widely studied within 
the Ginzburg-Landau framework [5J [TT] . They have the special structure 

w(x) = u(\x\)g ( xeR N 



where u is a solution of the following ordinary differential equation in 



pN 



N-l X K , 3 

-u = —u + u 



d|x| 2 |x| |x| 2 
u(0) = (46) 

and Xk is a characteristic constant. In what follows, we adapt Ginzburg-Landau techniques for (|46[) to the 
ordinary differential equation (|20p in the limit t — > oo to establish uniqueness and global monotonicity of h* 
in (j22)l . As will be demonstrated in Section [5l there are important technical differences between ([20]) and the 
Ginzburg-Landau formulation (|46[) and in general, Ginzburg-Landau results do not readily transfer to the 
Landau-de Gennes framework. In this sense, one could also refer to the t — > oo limit as the Ginzburg-Landau 
limit. 

Lemma 4.1. J^/ For all t > 1 and any solution Q of the Euler- Lagrange equations A15\) , we have the 
following global upper bound for the gradient - 

l|VQ|| L =o (fl(0)Ji) ) < C (47) 

where C > is independent oft. For the radial-hedgehog solution Q*(r) = *J\h*(r) (r (g> r — ^i), this 
implies that for t > 1, we have the following inequality 

|VQt=(^) 2 + ^<C 2 W€[0,R] (48) 

where C is again independent oft. 

Proof: The proof of Lemma 14.11 can be found in [5] where the authors show that a solution u of the 
elliptic system 

-Au = f onOc E n 

satisfies 

|Vw(r)| 2 < C\\f\\ L o* m \\u\\ L o* m Vrefi. 
In our case, we apply this result to the system (TT5|) . noting that Q* is a solution of (|15l) . 

/ = -Q y - - ^f 1 ( QifcQ ^ - T tr(Q2) ) + 2J r Q ^ (trQ2) 

for each i, j = 1 . . .3, ||Q*||l<x>(q) < 1 from the bounds in J24]) and ^ < | anc } V A < 31 f or t>l.U 
Lemma 4.2. In the limit t — > oo, we have that 

limr 2 ^=0. (49) 

r— >oo (XT 
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Proof: In the limit t — > oo, we solve the ordinary equation (|20[) on an unbounded domain i.e. the 
boundary conditions (p?Tj) become 

h(0) = h(r) ->lasr->oo. (50) 
From the bounds (f3"3")h we deduce that for r sufhciently large, 

a (3 

h*(r) = 1 + er(r) where < cr(r) < — - (51) 

for positive constants a,/3 independent of t. These bounds imply that h*(r) — » 1 uniformly as r — oo and 
from Proposition 15 . 21 in the next section, this implies that 

dh* da 
dr dr 

for r sufficiently large. 

We use PD]) to obtain an ordinary differential equation for 8 — as shown below : 



r 2 —-^ = -4r : 



d 2 <5 . dS 



dr 2 ' dr ' 4 <5 + 2r (fc* 8 - fc*) + r 2 (*h* 2 - l) 5 (52) 
where (5 — ► as r — ► oo. We can then use differential inequalities as in [3] to deduce that 

where 7>2/3>0isa positive constant and (3 has been in defined in (j5Tj) . Since /3 is independent of t and 
r, 7 can be chosen to be a positive constant independent of t and (I53[) implies that 

l imr 2^ =0 . (54) 
r— >oo ar 

Lemma l4~2l now follows. □ 

Proposition 4.1. TTie ordinary differential equation 

d 2 h 2 dh 6h , , ~ 3h + o 2 , 

+ - = -/i + /r 3 + — ± /i 3 -/i 2 (55) 

ar- 2 r dr r z t 

subject to the boundary conditions 

h(0) = h{r) -+ 1 r ^ oo (56) 
/las a unique non-negative solution in the limit t — > oo. 

Proof: Let /ii and be two different non-negative solutions of (|55[) subject to the boundary conditions 
dSSJ) i-e. 

f + - h ^-- 2 +{l-h\) + ^h l {l-h\)= Q 
h\ r hi r z x ' t x ' 

^ + -^-^ + ( 1 -^) + ^ 2 (l-^)=0 (57) 
«2 r h% T t ' 

where h x — h x — etc. We subtract the two equations to get 



= (l + (hi -h\) + ^ (hi - h 2 ) . (58) 



h± + !h _2 ( !h _!h\ _ fi , 3h + \ ^ ,2^ , 3h+ 

hi h 2 r \ hi h 2 
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Following the methods in [T], we multiply both sides of ([55} by r 2 (h\ — /i?,) and integrate from r = to 
r = R to find 

— n<2 \ ar + 

^ 2 „2 



fta)V(Ai + ha) dr 



/'I 



= - h2 J '* - 7-2/11 [t - hi ) "° ■ (59) 

Taking the limit R — > oo and using (14T?)) . we have that 

^£ r2 (fe^ ; ) 2 dr+ C r2 (^ h[ - k '>) 2 dr+ 

+ 1 R i 1 + ^ ^ ~ ^ 2 dr ~ I R {kl ~ k2)2 7-2 {kl + h *)dr = 0. (60) 

From (IBUl) . we deduce that 

£ (l + ^) r 2 - h 2 2 ) 2 dr-^± J" {hi - h 2 f r 2 (hi + h 2 ) dr -> R -> oo. (61) 

We first make the elementary observation that 3i?i e [0,iZ] such that 

hi(r) + h 2 (r) > 1 Vr > ifc. (62) 

The inequalities (1551) are true for any solution of (1551) subject to the boundary conditions (j5"6"j) [TT1 U]. 
Therefore, we have for t > 1, 

2 2 

</u(r),ft a (r)< _ 2 (63) 

^ < i? 2 < 14 (64) 

i.e. i?i can be bounded independently of t for t > 1. 

We partition the integral contribution in (j60|) into two sub-intervals [0, Ri] and respectively. 

/* f 1 + ^f) 7-2 ^ " ^ dr ~^rJ R {hl ~ 112)2 r2 {hl + h2) dr = 

= J" 1 (l + ^) r 2 (hf - h 2 ) 2 ^+r 2 (hi h 2 f (hi + h 2 ) dr + 

+ I* r 2 (h\~h 2 ) 2 + r 2 ^±(h 1 -h 2 f(hi + h 2 )[h 1 + h 2 -l] dr (65) 

JRi t 

and note from ([6"2")) that 

[ R r 2 (h 2 -h 2 ) 2 +r 2 ^±(h 1 -h 2 ) 2 {h 1 + h 2 )[h 1 +h 2 -l] dr> f ? r 2 (hf - h 2 f dr. 

JRi t JRx 



r 2 + 14 - w ' w - r 2 + 6/17 

and 
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Claim: For t sufficiently large, 



V { 1 + ^r) r2 ^ ~ hl)2 ^r r2 {hl ~ 712)2 {hl + M dr > \ 1 1 r2 ^ ~ hl)2 dr - (66) 



Recalling that R\ can be bounded independently of t, we note that 

' 3k+ r 2 (h x - h 2 ) 2 {h, + h 2 ) dr < 7l ^±i?f < 72 



where 71 and 72 are positive constants independent of t. Therefore, the claim in (|66[) is equivalent to 

St > ^ ^2 . (67) 



C^(h 2 -h 2 ) 2 dr 



for a positive constant 73 independent of t. 
We note that 



Rl 9 7? 3 



so that as t — > 00, we have two possibilities - (a) J Q Rl r 2 (h\ — til) 2 dr = 0(1) as t — > 00 and (b) 
J Q Rl r 2 (h\ — ft, 2 .) 2 dr — o(l) as t —> 00. In case (a), the condition (|57)l is clearly satisfied for t sufficiently 
large and the claim follows. 



For case (b), we have 



f Rl 2 

/ r 2 (hi -h 2 2 ) dr^Oast^ 00. (68) 
Jo 



From Lemma 14.11 and the global bound (j24|) , we obtain 

|V (hj - h 2 2 )\ < D (69) 
where D is a positive constant independent of t. Consider ro € [0,i?i] and let 

\(hj - h 2 ) (r )\ = a >0. 



Then from we have that 



1(^-^1) W|>f re 



a ao 



and therefore 

r 2 (ft 2 -ft 2 ) dr > / £y dr > 74a 5 



4 

where 74 is a positive constant independent of t. Combining the above with f|68|) . we have that ao — ► as 
t — > 00 and hence 

h 1 (r) = h 3 (r) re[0,i2i] (70) 

since the choice of ro is arbitrary and we are interested in non-negative solutions. 
From (j6Tj) and (j62j) . we have that 



fl 'l + r 2 (h\ h 2 2 ) 2 dr - ™± J* (hr - h 2 ) 2 r 2 (h x + h 2 ) dr > 

'' ' \ + ^t) r 2 (hf -til) 2 dr-^± fa - h 2 ) 2 r 2 (h, + h 2 ) dr + [* r 2 (h 2 - h 2 2 f dff.l) 
t J t Jo Jr.! 
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For case (a).(|6"6"]l holds and (fTTj) can be written as 

I rRi rR 
2 



r (K ~ K) dr + I r [K - K) dr < 

n 



< 



j (l + ^tt\ r 2 (hl-hl) 2 dr-^J^ (h, - h 2 ) 2 r 2 (h, + h 2 ) dr ^ t^oo (72) 



from which we deduce that 

(h\-hl) 2 = Q r€[0,R] (73) 

or equivalently 

hi{r) = h 2 (r) re[0,R} (74) 
For case (b), we have established in (170)) that hi(r) — h 2 (r) r 6 [0, Ri) and hence 

J* 1 + *!±±?j r 2 (hi - h 2 2 ) 2 dr-^± J*' - h 2 f r 2 [h x + h 2 ) dr -> t -> oo. 

From (fBTj) . we deduce that 

/ r 2 (h\-h 2 2 ) dr^O t^oo 

JR t 

and hence, 

h 1 {r) = h 2 {r) re[R u R] 
Combining the above with (|70p . we have that 



h 1 (r) = h 2 (r) re[0,R] (75) 

in case (b) too. Proposition 14. 1 1 now follows. □ 

We, next, illustrate the applicability of shooting arguments to the ordinary differential equation (|55[) in 
the limit t — > oo, so that the corresponding boundary conditions are (l5B|) [4|[5]. 

From Proposition 13. 11 we have that for any solution h of (|55|) subject to the boundary condition 

h(0) = 0, 

3 a constant a 2 such that 

/i(r) - a 2 r 2 r -)• 0. (76) 

Given 02, we denote the corresponding solution by h(a 2 ,r). We are interested in non-negative solutions and 
hence, we take a 2 > 0. By analogy with [4, 8 , we call a 2 the shooting parameter. We consider three different 
classes of solutions 



• V = \a 2 > 0; 3z e (0, R a ) such that dh ^ ,z) = oj 

• Q = {a 2 > 0; dh( 2" z) > and h(a 2 ,r) < 1 for all r > o} 

• Tl = {a 2 > 0; > OVr e (0, R a ) and max r£(0iJio) ft (a 2 , r) > l} 
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where R a is the maximal interval of existence of the solution h(a,2,r). Clearly 
and 

VUQUK = (0,oo). 

Our aim is to show that V and 1Z are non-empty and open. Then, Q is also non-empty. We have a 
unique solution of the ordinary differential equation ((55)) subject to the boundary conditions (l56l) in the 
limit t — > 00. Therefore, if we can show that a-i g Q implies that the corresponding h(a2,r) is a solution of 
(1551) and (|5l?|) . then we have that ft* e Q in the limit t — > 00 and hence, > for all r > i.e. we have 
global monotonicity in the limit t — > 00. 

It is evident that a solution of (|55D subject to the boundary conditions ([51))) cannot belong to 1Z owing to 
the global bounds (l24l) . It remains to rule out the possibility 82 £ P. We start with an elementary lemma. 

Lemma 4.3. If ' a-i G Q, then h(a2,r) is a solution of 155\) subject to the boundary conditions Ii56\) . 

Proof: The proof closely follows the methods in [3]. Since h(a,2,r) is monotonically increasing (from the 
definition of Q) and is bounded above by 1, b = linv^oo 6(02, r) exists and b g (0, 1]. Hence, to finish the 
proof, we need to show that 6=1. In fact, if b < 1, then as r —> 00, (|20| can be written as 

_1 ( =Qb + r : 



dr \ dr J 

so that 



^±(b 3 -b 2 )^j 



dh 6b r / , , 3h + , , 

contradicting the hypothesis that ^7 > for all r > 0. Therefore, 6=1 and Lemma |4~31 follows. □ 
Next, we need to show that Q is non-empty. For this, we need 

Lemma 4.4. The set V is not empty; more precisely, there exists a positive constant m such that (0, m) C V. 

Proof: The proof closely follows the methods in [4)[8]. For any 02 > 0, set 

h(a 2 ,r) 

w(a2,r) = ; (77) 

0-2 

then w satisfies the following ordinary differential equation from (|55t 

d 2 w 2 dw w 2 n 36.+ / 2 3 2\ 

—-5- -I : 6— + w = a 2 w H (a 2 w — a 2 w ) 

dr^ r dr r z t 

w(a 2 ,r)-r 2 r 0. (78) 

Then as 02 — > 0, w(a2, r) — > w(0, r) where w(0, r) is the solution of 

d 2 w 2 dw g w ^ 
dr 2 r G?r r 2 

u;(0,r)-r 2 r -> (79) 
and the general solution of this ordinary differential equation is 
C 

w(r) — — (— 3cosr — 3r sinr + r 2 cosr) + C 2 r 2 (— 3sinr + 3r cosr + r 2 sinr) (80) 

for arbitrary constants C\ and C2. From (|80|) . we deduce that w(a2,r) has oscillatory behavior as a 2 — >• 
and hence, so does ft(a2, r) = a 2 w(ci2,r). This completes the proof of the lemma. □ 
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Lemma 4.5. The set V is open. 



Proof: The proof of Lemma 14.51 closely follows the methods in [H [5] and we reproduce the proof for 
completeness and to illustrate the technical differences. 
For a 2 £ V, define 



Z o(a 2 )=inf|re(0, J R Q ); ^^1=0 

i.e. zo(a 2 ) is the smallest stationary point of h(a2,r). We can show that 

d 2 h(a 2 ,z (a 2 )) 



The definition of zo(a 2 ) implies that 

dh(a 2 ,z (a 2 )) 
dr 



dr 2 



= and 



< 0. 



d 2 h(a 2l z a (a 2 )) 
dr 2 



< 



since we are interested in non-negative solutions. 

From the governing ordinary differential equation (|55[) . we have that 



and note that 



If 



d 2 h 2 dh 
dr 2 r dr 

d_ 

dr 



h^ + h 2 -l+^{h 2 -h) y j<0atr = z {a 2 ) 



X + h 2 -l • ih+ 



^± (h 2 h) 



< at r = zo(a 2 ). 



then 



This implies that 



d 2 h(a 2 ,z Q (a 2 )) _ 
dr 2 

A + h 2 - 1 + ^± (h 2 - h) = at r = z (a 2 ). 
+ h 2 - 1 + ^± (/i 2 -/i)>0onr£ [z (a 2 ) - (5, z (a 2 )) 



for some 6 > 0. On the other hand. 

r 2 dr 



, d/i 



= h 



3h: 



(h 2 - h) 



from (EQl) so that 



r 2 dr 



, d/i 
dr 



> on r £ [z {a 2 ) - 5, z (a 2 )) . 



This in turn implies that 



z 2 (a 2 ) 



dh(a 2 , z {a 2 )) 



> {za{a 2 ) - 8) 



2 dh(a 2 ,z (a 2 ) - 5) 



> 



(81) 



(82) 



(83) 



(84) 



dr dr 

(since h r > for r £ (0, zo( fl 2 )) from the definition (|8ip). contradicting the definition of zo(a 2 ). Hence, ([8"2"j) 
holds. 

Finally, we note that for any oq £ V, by the Implicit Function Theorem and (|82|) . there exists a smooth 
function y(a 2 ) defined in a neighbourhood of ao such that y(ao) = zo(ao) and dh ( a2 ^y( a2 ^ = o. Hence, T 5 is 
open as required. □ 
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Lemma 4.6. The set 7Z is non-empty and open. 

Proof: The proof of Lemma [4T61 closely follows the methods in [5]. We introduce the function 

— 1/3 

v(r) = bh(a2,br) where b = a 2 ■ (85) 

Then one can check that v satisfies the following ordinary differential equation 

d 2 v 2 dv 6v / 3/i+\ o , 9 3/i+, 9 „ ,„„. 

+ f-[l + —±)v 3 + b 2 v + —±bv 2 = 86 

cfr^ r dr r z \ t J t 

with 

v(r) ~r 2 r -> 0. (87) 

If wc let 6 — > 0, then the limiting problem is 

d 2 v 2dv 6v ( 3h> 



dr 2 r dr r 2 \ t 
v(r) ~ r 2 r — > 0. 

From the hypothesis, we have that v, ^ > for 7- > 0. We claim that there does not exist I > such that 
liirv-^oo v(r) = I. We prove the claim by contradiction. Assume 31 > such that lim r _>oo v(r) = I. Then 
(1551) implies that 

*( 1 >to)„ a + 1 *( 1+ *!l± ]l * r^oo 
dr \ dr J \ t 



so that 

dv 61 r f 3/i+ 
dr r 3 \ t 



/ 3 r 



Therefore, w(r) >> I for r sufficiently large, which contradicts the hypothesis. The other possibility is I = 
but this contradicts the definition of 1Z which requires that ^ > for all r > 0. Therefore 

v(r) — > oo as r — > oo. (89) 

Consequently, h(a2,r) is large when ai is large enough and the set 1Z is non-empty. By the continuous 
dependence of h on 02 and the definition of 1Z, we deduce that 7Z is open. □ 

Lemma 4.7. The set Q is non-empty. 

Proof: This is immediate from Lemma 14.41 and 14.61 We omit the proof for brevity. □ 

Proposition 4.2. The function h* in A22\) is monotonically increasing in the limit t — > 00. 

Proof: From Lemmas I4.3I4.4I4T51 and 14.71 we have that there exists a a-i £ Q such that h(a,2,r) is a 
solution of (|5"5"|) subject to the boundary conditions (|5"6"|) . From Proposition 14. 11 we have that (|5"5"|) and (JSHJ) 
admit a unique solution h* in the limit t — ¥ 00. Hence, we deduce that the corresponding shooting parameter 
a* G Q i.e. h* is monotonically increasing everywhere away from the origin. An immediate consequence of 
this global monotonicity is < h*(r) < 1 for r € (0, 00). □ 
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5 The L ->> limit 

Consider the Landau-de Gennes energy functional in ([3]) 

l£0[Q]= I §|VQ| 2 + / B (Q) dV 
Jb(o,r) z 

in the limit L — > + . Let Q L be an arbitrary solution of the corresponding Eulcr-Lagrange equations 



df B 1 df_ 



(where | ecitk ^ 3 * s a Lagrange multiplier accounting for tracelessness) that satisfies an energy bound of the 
form 

leg [Q L ] < LC(a 2 , b 2 , c 2 , n) VL > (90) 

where C does not depend on L in the limit i — >• + . Examples of such solutions include global energy 
minimizers and the radial- hedgehog solution (see (j2U)l ). Then it is intuitively clear that as L —> + , such 
solutions will be almost like the bulk energy minimizers in ([5]) and (jSJ) (since the energy bound implies 
Xb(O-R) /b(Q) ^ — LC (a 2 , b 2 , c 2 , ft) as L -> + ), with the elastic energy density being dominant in the 
vicnity of defects and interfaces. This has been rigorously established for global energy minimizers in [14] 
and for solutions satisfying the energy bound in (f9"0")l . in |16| . 

Consider the ordinary differential equation for h* in ((20]) and the boundary conditions ((21"]) in the limit 
L -> 0+ i.e. 

d 2 /i 2 dh 6h ^ 3 3/i+ / 3 2 \ 

dr 2 r dr r 2 t 

h(0) = h{r) -> 1 as r -> oo (91) 

since the ball radius is inversely proportional to the correlation length £ defined in ([9]). There is an important 
difference between (|9Tj) and ((55)) which focusses on the limit t — > oo. In ((9"Tj) . the term (/i 3 — /i 2 ) is not 
necessarily much smaller than the term —h + h 3 on the right-hand side of (|9~T1) . Hence, ((9"Tj) does not have 
the Ginzburg-Landau structure as in (|46l) and we do not have analogous uniqueness and global monotonicity 
results for the radial-hedgehog solution in the L —> limit. We can, however, utilise Ginzburg-Landau 
techniques to understand the far-field properties, away from the isotropic defect core, in the L — > + limit. 

We first recall an important result from (14) and |16] regarding the uniform convergence of h* everywhere 
away from the isotropic core. 

Proposition 5.1. \1J\ \16§ For fixed t > 1 and L sufficiently small, there exists R > such that h*(r) — > 1 
uniformly for all r > R. 

Proposition 5.2. Let h* be a global minimizer of I in \18\) , in the class Ah, in the limit L — > 0. Then there 
exists Ri > such that h* is monotonically increasing for all r > R\. 

Proof: From Proposition 15. 1[ we have that there exists Rq > such that h(r) > h for r > Rq. Consider 
the right-hand side of (|9"T|) and define 

F(h) = h 2 -l + ^±(h 2 -h). (92) 
Then F(l) = and F' (h) > for h > \. 
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We prove Proposition 15.21 by contradiction. We assume that there exists ro > R, where R is defined in 
Proposition 15. 11 such that 

— \r=r — U. 

dr 



There are three possibilities for 44-| r=nj i.e. (a)^4-| r: 
Consider case (a). Then we have from (l9"TT) that 



-r = 0, (b) ^\ r =r < 0, and (c) ^| r=ro > 0. 



d 2 h* 2 dh* 



r dr 



dr 2 



= h* 



= at r = ro- 



(93) 



Secondly, 



from which we deduce that 



d_ 

dr 



6 



< at r — rg, 



F ( h *) + ^ >0 re(r -S,r ) 



for some S > 0. We deduce from (1911) that 



(94) 



d 
dr 



,dh* 
dr 



> r € (r - 6, r ) 



so that 



|r=r > (^ ~ S) 2 —\ ro -g 



Since ^r\ r=rg = 0, we deduce that ^fe-| ro -5 < 0. This necessarily means that there exists a local minimum 
at r = r\ > ro, since < h*(r) < 1 and /i* — »• 1 as r — > oo. We, therefore, have 

d 2 h* , 



dr 2 



or equivalcntly 



> 



6 



F(h*(n)) + -j > o 



But 



F(^(n)) + -! < ^*M) + -! = 



since F'(h) > for /i > \ and /i*(n) < h*{r ). This gives a contradiction and we deduce that 



dr 2 |r=r 



=r ¥= 0. 



Case (b): We assume that ^r\ r =r < and h*(ro) > \ i.e. we have a local maximum at r = ro. The 
local maximum must be followed by a local minimum at r — r± > ro, since < h*(r) < 1 Vr > and 
h* — > 1 as r — > oo. Thus, 

F(ft*(n)) + 4 > o 

by definition of a local minimum from (1931) . However h(r±) < h(ro) and 

F(h*(n)) + ^<F(h*(r )) + ^<0 



yielding a contradiction. 
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Case (c): We assume that 



drh* i 



> 0. Then > for r > r > R where R has been defined 



in Proposition 15.11 since the previous arguments show that we cannot have a point of inflection or a local 
maximum for r > R. Then we set R\ in Proposition 15 . 21 to be R\ = r$. Proposition 15.21 now follows. □ 

We use the far-field monotonicity established in Proposition 15.21 to derive an explicit far- field expansion 
for h* as r — > oo. This expansion is valid in both the L — > and t — > oo limits. 

Proposition 5.3. Let h* be a minimizer of I in \18\) in the space Ah, for a fixed t > 1, in the limit L — > 0. 
Then h* is a non-negative solution of the following singular ordinary differential equation 



d 2 h 2 dh 6h 



3h, 



dr 2 r dr r 
subject to the boundary conditions 

h(0) = and h(r) ^ 1 as r 
We have the following far-field estimates 

6-r 2 h*(l- h*) (1+ [1 + 



oo. 



d 2 h* 




dh* 




dr 2 


+ r 




+ 


dr 





3h, 



o(l) r 



oo. 



(95) 



(96) 



(97) 



Proof: The proof of Proposition 15.31 follows some of the methods described in a recent paper [TS] on 
Ginzburg-Landau theory for three-dimensional domains. 

The bounds (|33[) are valid in the L — > limit. In fact, they are valid in any limit which translates to an 
unbounded domain in terms of the dimensionless variables in ([§]). In particular, they imply that 

a 
l_«<*. (r) <l_£ 

as r — > oo, as shown in ((5"Tj) . for positive constants a, independent of L. As demonstrated in (|49|) . this 
implies 



2 dh* 
lim r — — 

> — >oo dr 







and hence 



r dh * n 
hm r— — = 0. 

r^oo dr 



(98) 
(99) 



For any k € (0,1) fixed, we multiply (|9"S"|) by r 2 , average over (fci?*, i?*), take the limit i?* — > oo and 
obtain 



1 



(1 - fc)i?* dr 



,dh* 
dr 



dr- 



(1 - k)R* J kR 



R, 



r~li 1 ( /■') ( ] - /rf/'i) | I + ( 1+ '^j- ) h j ,//• = 



G 



R, 



(1 - k)R, 
(100) 

In the limit L — > + , h* — > 1 uniformly as r — >■ oo from Proposition 15.11 and using (|98[) . we obtain the 
following sequence of inequalities 



/i*(r) dr. 



HmsupA; 2 ^(l-/i*(J?*))) ( 1+ ( 1 + — ) /i*(E*) ) < 6 < hminf Rl (1 - h*{kR*)) ( 1 



It immediately follows that 



(l-r(r)))(l+(l + ^± )^(r) 



(101) 
(102) 
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uniformly in the limit r — > oo. 

Finally, using the estimates and (|102|) in ([95)1 , we deduce that 



r 2 



c?r 2 



->■ (103) 



uniformly in the limit r — > oo. Proposition [531 now follows. □ 
One immediate consequence of ((57)) is that 



h*(r) = l y-^ r^oo. (104) 



(2+^) 



Although this information is qualitatively contained in (l5~Tj) . (| 104)) is a stronger result since it is an exact 
expression that captures the effects of geometry and the temperature on the far-field structure. Further, 
(|104[) yields estimates for the higher-order derivatives of h* as r — > oo and this information cannot be 
immediately inferred from (|51[) . 

Proposition 5.4. Let Q L denote the radial-hedgehog solution in V2S\) for a fixed L > and fixed t > 1. 
Tfcera {Q L } -> Q° W 1 ' 2 (B(0,R);S) as L -> 0+, w/iere Q° = t/| - |l) f^F- 4si4 0+ the 
radial-hedgehog solution Q L is stable against all small far-field perturbations P i/iai satisfy 

P(r)=0 reB(0,i? L ), 

P(r)=0 redB(0,R) (105) 

and i?L is sufficiently large. In other words, Q L is locally stable against perturbations which are localized 
outside the isotropic core around the origin, in the limit L — > + . 



Proof: Consider the dimensionless free energy in (flO| 

HQ] = I ^|VQ| 2 + / B (Q) dV 



where 



f b = -itrQ 2 - ^±trQ 3 + |- (trQ 2 ) 2 + C(t) > 



and 



/fl(Q) = 0^Q = y- ln®n-- 
We consider a small perturbation P that satisfies (|105|) . Define 



We compute the second variation 



A direct computation shows that 



Q^=Q^+ eP y 0<e«l. (106) 



d 2 I[Q<] , 



^[Ql, _ /■ |vp| 2 , Pf, 

df 2 U=0 - / |VP| + o ni o n L 

« e Jb(0,R)\B(0,R l ) U^iiV^m 



e=o = / |VPr + art£ 1" 7 PijPp, ^ (107) 
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since P = on B(0,R L ). 

From Proposition 15.31 we have that as L — > 0, 



Q L (r) 



1 



r r I 

-®-~ o 
r r o 



r — > oo. 



Equation (|108[1 can be written as 

Q L (r) = Q°(r)- ' 



r r I\ . 1 . 



where Q° is a bulk energy minimizer by dchnition (see preceding comments). 



We perform a Taylor expansion of 



Q 2 .fB 



around Q° to obtain 



d 2 fs 

Finally, we note that 



d 2 f B 



P P _ 7. ( r ^ 



'« ^pq 



d 3 f E 



3 J dQ%dQ%dQ% 



d 2 h 



aQ°.9Qo 9 PyPp9>7?(t) 

where r\ is independent of r (since Q° is a global minimizer of fs by definition) and 



_7 / r a r /3 



3 ; aq^SQ^Qo, 



^ r> - RL 



(108) 



(109) 



o r — > oo. 



(110) 



since the derivatives of fs can be bounded independently of Rl- Combining the above, we have 



d 2 fB n n 



> o 



for a fixed t > 1 and Rl sufficiently large. Substituting (jlllj) into (|107p . we deduce that 

d 2 I[Q € 



de 2 



-U=o > 



for perturbations P satisfying (|105[) . The positivity of the second variation ensures that the radial- hed 
solution is locally stable against perturbations P satisfying (|105[) . in the limit L — >• + . Proposition 
follows. □ 



(111) 



(112) 

;ehog 
now 



6 A general stability result 

In the previous sections, we have demonstrated that the radial-hedgehog solution is unstable with respect to 
biaxial perturbations, localised around the isotropic core, for R and t sufficiently large. As already stated, 
the R — > oo limit is equivalent to either the L — > or t — > oo limits or both. We have also shown that the 
radial- hedgehog solution is locally stable with respect to far-field perturbations in the R — > oo limit. We 
conclude by deriving a general local stability result that is not restricted to the limits t — > oo or L — > 0. 
We note that Proposition 16. II is known from numerical investigations (see [241 ?, ?]) and we present a proof 
partly for completeness and partly this proof gives greater insight into how the elastic constant, temperature 
and ball radius collectively affect stability properties. 
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Proposition 6.1. Let B(0,R) denote a ball of radius R centered at the origin in R 3 . The corresponding 
radial-hedgehog solution is stable against all small, smooth perturbations of the form 

Q = Q* R + eP (113) 

where e € K, |e| << 1, P G S and P = on dB(0, R), provided that the radius R is sufficiently small i.e. 

R 2 <\ \ K— I • (114) 



1 



+ 



In terms of the original variables defined in {P|) 7 Ijll4\l is equivalent to 



^<M77^I (115) 



where £ is the correlation length defined in 0). 



Proof: The results in Proposition 12.11 are true for any R > i.e. for every R > 0, we are guaranteed the 
existence of a radial- hedgehog solution of the form (f2"2")l . that satisfies the energy bound (|2"3"|) and the 
inequalities Consider the dimensionless free energy in (TT^)) and introduce the change of variable 



t 

r = — 
R 

so that the free energy becomes 

I[Q] = f V r t Ul V Q| 2 - ? tr Q 2 - ^^^trQ 3 + ^±i? 2 (trQ) 2 + R 2 C(t)\ f 2 sin Odfdddcj,. 
Jo Jo Jo [ 2 2 i 2t J 

(116) 

We consider small perturbations 

Q^Q^ + eP 0<e«l (117) 
such that P = on dB(0, R). We compute the second variation of the Landau-de Gennes energy functional 

^/[Qe]|e=0 = fj [ | |VP| 2 - i? 2 |P| 2 - + ^ [8 (Q* R ■ P) 2 + 4|P| 2 |Q^| 

(118)' 

where dV = f 2 sin 8dfd9dcf). 

We, next, make an elementary observation 

PyP JP Qfl pl = h*(r) [v t P^r p P p3 /r 2 - |P| 2 /s] < ||P| 2 



so that 



^/[qji„.> r r r w-^ipp-^ituvippi-a.*** <ii9) 



" c JO Jo Jo 

We note that 

IPI 2 > ( — 
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and use the following inequality from [SJ [T2] 

for a real- valued function a defined on the interval [0, 1]. Substituting the above inequality in (|120l) . we have 
that 

|^[Qe]|e=0 > fj J* | \\P\ 2 |P| 2 i? 2 + ^j^j | sin0 df d0 d^ (120) 
since r < 1. It follows that 

if 

or equivalently if 



^[Qe]|e=0>0 

R 2 < ~ ^— (121) 



y^i (i22) 



1 + 



+ 



t / 

where i? rea / = t/j-R from ([§]) and (fTTj) . Proposition 16.11 now follows. □ 



7 Discussion 



This paper aims to build a self-contained and rigorous mathematical framework for the study of the radial- 
hedgehog solution within the Landau-de Gennes theory for nematic liquid crystals and elucidate the analogies 
between the mathematical formulation of defects in the Landau-de Gennes framework and defects in the 
Ginzburg-Landau theory of superconductivity. These analogies need to be highlighted in the applied math- 
ematics literature, so that mathematical techniques from other branches of condensed matter science can be 
effectively used in the context of liquid crystals. We study radial-hedgehog solutions on spherical droplets 
subject to homeotropic anchoring or strong radial anchoring conditions and define a radial-hedgehog solution 
to be an energy minimizer within the class of spherically symmetric uniaxial solutions as demonstrated in 
Proposition 12.11 We consider two different limits in this paper: (a) the low-temperature limit t — > oo and 
(b) the vanishing core limit L — > + . For completeness, we summarize the validity of the different results 
in this paper in different parameter regimes. Proposition 12. II and Proposition 13. II are valid in all parameter 
regimes i.e. we are always guaranteed the existence of a radial-hedgehog solution that satisfies the energy 
bound (|2U)) and whose scalar order parameter is constrained by the inequalities (|24l) . We have a single 
isolated isotropic point at the origin by definition and we always have a series expansion near the origin that 
involves only even powers of r as r —¥ 0. However, the bounds (|33[) are only valid in the limit R — ¥ oo, where 
R is the re-scaled ball radius. Recalling the definition of the dimensionless variables in (O and (TTT|) . the 
limit R — > oo is equivalent to either the low-temperature limit or the vanishing core limit. Propositions 15.31 
and Proposition 15.41 are valid in the L —> limit, whereby it is difficult to establish rigorous results about 
the defect core but the governing equation has a Ginzburg-Landau structure away from the origin. This 
Ginzburg-Landau structure gives us a good grip on the far-field properties i.e. uniform convergence of the 
scalar order parameter away from the origin, far-field monotonicity and explicit far-field expansions for the 
scalar order parameter. The limit t — > oo has a Ginzburg-Landau structure (see (pES"]) and (1461) ) and we can 
exploit Ginzburg-Landau techniques to prove global properties. Propositions 14.11 [4.21 15.31 15.41 and 13.21 hold 
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in the t — > oo limit and we demonstrate the manifestation of biaxial instabilities localised near the isotropic 
core in this regime. Proposition 14. II is an example of how Ginzburg-Landau techniques can be used to prove 
results on multiplicity of solutions and in Proposition ^. 21 we appeal to shooting arguments which have not 
been used previously in the Landau-de Gennes context. Proposition 16. II is a general result that identifies a 
relationship between the elastic constant L, the reduced temperature t and the ball radius R that guarantees 
local stability of the radial-hedgehog solution against all perturbations. 

In [3J, H. Brezis postulated the following problem in the context of Ginzburg-Landau theory for super- 
conductors: for maps u : R 3 — > R 3 , is any solution of the system 



satisfying |u(r)| — > 1 as |r| + oo (possibly with a good rate of convergence) and deg^u — 1 of the form 



for a unique function / vanishing at zero and increasing to one at infinity. In j!8j . the authors show that 
every non-constant local minimizer of the Ginzburg-Landau energy functional associated with (|f 231) , 



is of the form (|124[) . up to a translation on the domain and an orthogonal transformation on the image. For 
nematic liquid crystals, the corresponding problem translates to: is any uniaxial solution of (fT5)l necessarily 
of the form (|22]) i.e. are radial-hedgehog solutions the only possible uniaxial solutions of the system (fl~5)) in 
R 3 ? If so, then we will have a complete characterization of all admissible uniaxial solutions and the interplay 
between uniaxiality and biaxiality can be partially understood in terms of the comparitively tractable radial- 
hedgehog problem. We expect that the methods in [Tj5] will not readily transfer to the Landau-de Gennes 
framework and there will be analogies only in certain parameter regimes, such as the t —¥ oo limit studied 
in this paper. 

Finally, we compare our results with previous work in this area. In |24j . the authors carry out detailed 
numerical investigations of equilibrium configurations within spherical droplets subject to strong radial 
anchoring conditions and find that the radial-hedgehog solution only occurs either in very small droplets 
or very close to the nematic-isotropic transition temperature; the symmetry-breaking biaxial torus solution 
is energetically preferable everywhere else. This is consistent with Proposition 16.11 where we demonstrate 
local stability of the radial-hedgehog solutions for droplets with radius comparable to the nematic correlation 
length £. This is also consistent with Proposition ^. 21 where we demonstrate that the radial-hedgehog solution 
cannot be a global energy minimizer for large droplets in the low-temperature limit. In [22] . the authors 
work within the Lyuksyutov constraint, which requires that trQ 2 (r) = |s^_ for r G B(0,R), where s+ has 
been defined in ((BJ. They demonstrate that the radial- hedhegog solution is always locally unstable within 
the one-constant approximation for the elastic energy density i.e. when the elastic energy density is simply 
taken to |VQ| 2 , as has been done in this paper. This resuly is evidently in agreement with Proposition 13.21 
and does not contradict Proposition 16.11 where we demonstrate local stability in balls of sufficiently small 
radius. The Lyuksyutov constraint is valid in the R — > oo limit or for balls of sufficiently large radius and 
hence, Proposition 16.11 is outside the remit of this instability result. The analogies of this work with the 
results reported in [TTl [9] have already been mentioned. The authors numerically study the stability of the 
radial-hedgehog solution as a function of the ball radius, reduced temperature and elastic constants in [51 IT?] 
( along with analysis of the biaxial instabilites as has been mentioned in Section [3]) and in the one-constant 
case, our results are qualitatively in agreement with the phase diagrams in [HI HZ]- It would be interesting 
to see if (|122p can yield a qualitative fit to the region of local stability obtained in [5J [T7J . While careful 




(123) 



U(r) = ~/(r) 



(124) 
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attention is paid to the effect of elastic constants in some of the previous work, we focus on the one-constant 
case. This is primarily because the one-constant case has a much more tractable mathematical structure 
than the unequal constant case and is the best paradigm for illustrating the efficacy of Ginzburg-Landau 
techniques in the Landau-de Gennes framework. The unequal elastic constant case will be considered in 
future work. 
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